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The collision of two atoms is an intrinsic multi-channel (MC) problem as becomes especially 
obvious in the presence of Feshbach resonances. Due to its complexity, however, single-channel (SC) 
approximations, which reproduce the long-range behavior of the open channel, are often applied 
in calculations. In this work the complete MC problem is solved numerically for the magnetic 
Feshbach resonances (MFRs) in collisions between generic ultracold ^Li and *^Rb atoms in the 
ground state and in the presence of a static magnetic field B. The obtained MC solutions are used 
to test various existing as well as presently developed SC approaches. It was found that many 
aspects even at short internuclear distances are qualitatively well reflected. This can be used to 
investigate molecular processes in the presence of an external trap or in many-body systems that 
can be feasibly treated only within the framework of the SC approximation. The applicability of 
various SC approximations is tested for a transition to the absolute vibrational ground state around 
an MFR. The conformance of the SC approaches is explained by the two-channel approximation for 
the MFR. 



I. INTRODUCTION 

The tunability of the interparticle interaction on the 
basis of Feshbach resonances, especially magnetic ones 
(MFRs) , marked a very important corner-stone in the re- 
search area of ultracold atomic gases. At ultracold ener- 
gies s-wave scattering dominates the atom-atom interac- 
tion, such that for large internuclear distances the elastic 
scattering properties are solely described by the s-wave 
scattering length Ogc [J. Its sign determines the type of 
interaction (repulsive or attractive) and its absolute value 
the interaction strength. In the presence of an MFR this 
parameter can be tuned at will by applying an external 
magnetic field. A wide range of experiments using MFR 
techniques has been carried out including the formation 
of cold, even Bose-Einstein condensed molecules [2], 0, Q 
or the realization of a Mott insulator phase with atoms 
in an optical lattice (OL) 

Experiments with ultracold gases are usually per- 
formed in external trapping potentials and over an en- 
semble of many particles. For tight trapping conditions 
the influence of the additional potential can become es- 
sential. For example, processes of molecule formation via 
photoassociation (PA) where two ultracold atoms absorb 
a photon and form a bound excited molecule [1, 0] can 
be more efficient, if performed under tight trapping con- 
ditions as they are accessible in OLs [Jl^, ■ 

However, the presence of a trapping potential or, 
worse, the existence of many-body effects is a great chal- 
lenge for the full theoretical description of an MFR, since 
all accessible spin configurations of the colliding atoms 
must be included, leading to a multi-channel (MC) prob- 
lem. For the case of s-wave scattering of two free atoms 
the separation in relative and center-of-mass motions, 
the formulation in spherical coordinates, and the con- 
tinuous energy spectrum make the numerical solution 
manageable. This changes, unfortunately, if an exter- 
nal potential couples the six spatial coordinates of the 
two colliding atoms and induces the need to find dis- 



crete eigenenergies This can be the case for atoms 
loaded in a cubic OL formed with the aid of standing 
light waves Furthermore, the theoretical mi- 

croscopic investigation of ultracold many-body systems 
is feasible only within the framework of the SC approxi- 
mation. Nevertheless, a good knowledge of two-body MC 
collisions should help in understanding the consequences 
of SC approximations which must be done when many- 
body systems are considered. 

Single-channel (SC) approximations allowed to study 
the influence of the scattering length as it results from an 
MFR for three-body collisions [la, [3 and in the presence 
of an external trap [^, [l^, However, to our knowl- 
edge it is not yet well established to what extent SC ap- 
proximations describe correctly the behavior of a coupled 
MC system, if more than one channel contributes signifi- 
cantly. The successful usage of (SC) pseudo-potentials to 
model, e. g., the atom-atom interaction in OLs [l^ shows 
that physical properties depending on the long-range be- 
havior of the open-channel scattering wave function, i. e., 
the scattering length age, are well described within the SC 
framework. For shorter interatomic distances in the or- 
der of the van der Waals length scale j3Q (/Jg = {^^-CqY^^ 
where /i is the reduced mass and Cg is the van der Waals 
coefficient) this is not necessarily the case. Here, all cou- 
pled channels contribute to the full wave function and 
affect processes, such as transitions to molecular bound 
states. For these distances, SC approximations cannot 
cover all details of the MC solution. As will be shown, 
some important aspects are, nevertheless, reflected and 
can be used to study processes of molecule formation in 
the presence of an MFR where MC calculations may be 
too laborious. A very systematic investigation of both 
short-range and long-range parts of the MC solutions 
against various SC ones is considered that is done in this 
work. 

The formation of ultracold molecules especially in 
deeply bound levels is currently of large interest. In order 
to associate them, the starting point is often a sample of 
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Feshbach molecules, obtained from ultracold atoms via a 
sweep of the magnetic field around an MFR 0, [l^, . 
These molecules are usually formed in high lying vi- 
brational ground states. Molecules in lower vibrational 
states and eventually in the absolute vibrational ground 
state are, however, favorable since they are more stable 
against inelastic collisions. The most successful scheme 
to access those molecules is the two-color stimulated Ra- 
man adiabatic passage (STIRAP) [2l| when the passage 
is realized using PA via an intermediate excited state. 
The dump photoassociation (DPA) process during which 
two ultracold atoms absorb a photon and form directly a 
ground molecule is in principle possible for hetcronuclear 
systems, although the yield is very small. 

It has been shown theoretically and for some cases even 
experimentally, that the PA and DPA yields can be sig- 
nificantly increased in the presence of an MFR [3, [l^, [23, 
SHE, [ill. For example in it was found that an 
SC scheme based on mass variation predicts the same en- 
hancement of the PA rate for almost all final states except 
the very high-lying ones and the ones at the PA window 
(for Use > 0). The reason for the enhancement was the 
increase of the absolute value for the initial-state wave 
function that occurs for large absolute values of age. As a 
consequence the corresponding Franck-Condon (FC) fac- 
tors and PA yields increase with |asc| (see Sec. Ill, G of [l^ 
for details). Noteworthy, a strong enhancement of the PA 
rate by at least two orders of magnitude while scanning 
over an MFR was predicted on the basis of a MC calcu- 
lation for a specific ^^Rb resonance already in The 
explanation for the enhancement given in [2^ is, how- 
ever, based on an increased admixture of a bound-state 
contribution to the initial continuum state in the vicin- 
ity of the resonance. This is evidently different from the 
reason for the enhancement due to large values of |asc| 
discussed in [l^]. This suggests that both seemingly dif- 
ferent systems appear exhibit a strong correspondence. 
One of the motivations of the present work is to clarify 
this observation. 

To mimic certain aspects of the MC wave function for 
studying molecular processes, SC approaches make use 
of a controlled tuning of system parameters such as the 
reduced mass [l^, van der Waals coefficients [l^l, inner 
wall [ll| of the interaction potential or the interaction 
potential in the intermediate range as is proposed in this 
work. Long-range scattering properties like the s-wave 
scattering length can be sensitive to even small changes 
of those parameters. To date, the justification of these 
systematical variations is mainly given by the broad va- 
riety of atomic species and their isotopes, each with dif- 
ferent parameter values. In this work it will be shown 
that by these variational approaches one is also able to 
reproduce changes of both long and short range collision 
properties of a given scattering system as it is induced by 
an external magnetic field in the proximity of an MFR. 

The general validity of the SC methods will be based on 
a two-channel (TC) approximation of the MFR [13, [ll]. 
This approximation is widely used to describe the phe- 



nomenon of an MFR and has been adopted to study 
many-body interactions [2^ and two-atom interaction in 
a time-dependent magnetic field S^Slj and in a struc- 
tured continuum induced by an OL [32|. The TC approx- 
imation reproduces many aspects of the coupled MC sys- 
tem. It allows to describe the complex PA transition pro- 
cess by just two free parameters, the maximal transition 
rate and the position of the minimal transition rate [33j . 
An analysis of the TC approximation reveals why SC ap- 
proaches can show an astonishing conformance with the 
coupled MC predictions. 

In order to compare concrete MC and SC solutions 
the exemplary case of ^Li and ^^Rb scattering is consid- 
ered and the relative motion of this system in a static 
magnetic field B is fully solved employing the i?-matrix 
method (34| . This system is of great importance by it- 
self for its large static dipole moment, which makes it 
interesting for applications in quantum information pro- 



cessing [35. 136i or the exploration of lattices of dipolar 
molecules 37l|. The applicability of the different SC ap- 
proaches is studied by considering the process of molecule 
formation by a direct PA of ^Li-^^Rb to the absolute vi- 
brational ground state in the presence of an MFR. We 
describe this process by using the exact MC solution and 
compare to different SC approximations. 

The paper is organized in the following way. In Sec. [Til 
a theoretical description of ^Li and *^Rb scattering is 
given and the TC approximation is briefly introduced. 
The possibility of SC approaches is motivated by con- 
sidering the results of a full MC calculation for differ- 
ent resonant and off-resonant magnetic field values. In 
Sec. mil diverse SC approaches arc introduced and their 
wave functions are compared to those of the full MC cal- 
culation. The direct dumping to the absolute vibrational 
ground state is considered in Sec. |TVl The prediction 
of the TC approximation is presented and MC and SC 
results arc compared. Finally, a conclusion is given in 
Sec. [V] All equations in this paper are given in atomic 
units unless otherwise specified. 



II. MULTI-CHANNEL APPROACH 

A. Hamiltonian 

The Hamiltonian of relative motion for two colliding 
ground-state alkali atoms - in the present case ^Li (atom 
1) and *^'^Rb (atom 2) - is given as [38[ 



H = t^+^(Vf +Vj^)+Vint 



(1) 



where T^ is the kinetic energy and ^ is the reduced mass. 
The hypcrfinc operator V 



hf 



ij and the Zceman 



operator = (7eS*j — jnij) • B in the presence of a mag- 
netic field B depend on the electronic spin s,- and nuclear 
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spin ij of atom j = 1,2. For the present system the val- 
ues of the hyperfine constants a^£, a^f, and those of the 
nuclear and electronic gyromagnetic factors 7„ and 7e 
are adopted from [s^. In Eq. ^ the central interaction 
Vint (R) between the atoms is a combination of electronic 
singlet and triplet contributions 



YUR)^VoiR)Po + Vi{R)Pi 



(2) 



where Pq and Pi project on the singlet and triplet compo- 
nents of the scattering wave function, respectively. The 
potential curve Vq (Vi) for the singlet (triplet) states 
of ^Li-'^^Rb in Born-Oppenheimer (BO) approximation 
were obtained using data from [40, [4l| and references 
therein. In [40| refined potential parameters such as the 
van der Waals and exchange coefficients, which we use in 
the following, were determined by a comparison of MC 
calculations with experimentally observed resonances. It 
is important to note that the MC approach considered in 
the present work is formulated in relative motion coordi- 
nates. This is based on the assumption that the centcr-of- 
mass and relative motion of two atoms may be decoupled 
and effects due to coupling may be neglected. Further- 
more, calculations of the pres ent work assume the BO 
approximation to be valid |42| . 

For the interactions present in Hamiltonian ([1]) the pro- 
jection Mp of the total spin angular momentum F = 
fi + /2 on the magnetic field axis is conserved during the 
collision. Here, fj = Sj + ij is the total spin of atom j. 
For a given Mp of the colliding atoms only spin-states 
with the same total projection of the angular momen- 
tum can be excited during the collision. If {|a)}Q, is a 
complete basis of the electron and nuclear spins of the 
Mi?-subspace, one may use the function 



MR) , 

R ' 



(3) 



in order to find the s-wave scattering solution of the sta- 
tionary Schrodinger equation with Hamiltonian ([1]) . This 
ansatz yields a system of coupled second-order differen- 
tial equations 



2^aR2 



Vc^iR) + E^{B) - E jMR) 

+ Y,Wo.'c.{R)M{R)^Q 



(4) 



where the channel threshold energies Ea , the channel po- 
tentials Va{R), and the coupling potentials Wc,'a{R) de- 
pend on the chosen spin basis and will be specified below. 

Depending on the spin basis, the scaled channel func- 
tions ipa{R) will be used in the analysis instead of the 
full channel functions ipa{R)/R, while the name "chan- 
nel function" is kept for convenience. 



1. Atomic basis 

If the two atoms are far apart from each other, the cen- 
tral interaction Vint(-R) may be neglected and the two- 
body system is described by the spin eigenstates | /j , m f. ) 
of each atom. In this atomic basis (AB) the collision 
channels \a) are written as a direct product of the atomic 
states \x) = l/i, m/JI/2, m/j). In this case the thresh- 
old energy E^{B) of channel |x) is given as the sum of 
Zeeman and hyperfine energies of the two atoms. The 
channel potential V^{R) in the AB is identical for all 
channels, 



V^{R) = V+{R) 



Vo{R) + Vi{R) 



(5) 



The long-range asymptote of V+ is described by an at- 
tractive van der Waals interaction, that in the present 
case of ^Li and ^^Rb atoms in their ground states is given 
as 



C 

n—3 



(6) 



with Cg 2543 a.u., Cg 228250 a.u., and Cm = 
25 645 000 a.u. The coupling between the channels in the 
AB is given as l^x'x(^) = (x'|Po - Pi|x)^-(^) where 

= (7) 

The exchange interaction Vex is in the long-range regime 
very well represented in the Smirnov and Chibisov 
form El 



VcKiR;Ja,a) = JqR'^ 



(8) 



In Eq. ([5]) Jo = 0.0125 is a normalization constant and 
a = 1.184 depends on the ionization energies of each 
atom. For a given magnetic field B the channel thresh- 
old energies E^ and coupling matrix W^^i are fixed and 
V^{R) describes how strongly the different channels \x) 
are coupled. 

The total energy E available to the system is the ki- 
netic energy, i. e., the energy at a time prior to the inter- 
action when particles are far apart from each other. Since 
the coupling vanishes exponentially, the channels in the 
AB are asymptotically uncoupled. If the threshold en- 
ergy of a channel either lies above or equals the total en- 
ergy available to the system, E^{B) > E, the channel is 
considered to be "open" , otherwise it is "closed" . With- 
out loss of generality we consider in the following an elas- 
tic collision where only the channel |ai) = |l/2, 1/2) |1, 1) 
with the lowest threshold energy is open. The thresh- 
old energy Ea-^ marks the zero point of the energy scale 
throughout the paper. 



2. Molecular basis 

Another possible choice of the spin basis of the 
channels \a) is the molecular basis (MB) |^) = 
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TABLE I: Atomic and molecular basis states of the ^Li-*^Rb 
system for the manifold of states with Mf ~ 3/2. 



index \x) atomic basis index |^) molecular basis 



l«l> 


11/2,1/2)11,1) 


\Si) 


|0,0)|l,l/2) 


l«2> 


13/2,1/2)11,1) 


\S2) 


|0, 0)10,3/2) 


l«3> 


13/2,3/2)11,0) 


\Ti) 


|l,-l>|l,3/2) 


|a4> 


11/2,1/2)12,1) 


\T2) 


|1, 0)10,3/2) 


las) 


11/2,-1/2)12,2) 


ITs) 


|l,0)|l,l/2) 


lae) 


13/2,3/2)12,0) 


\T4) 


11,1)1-1,3/2) 


l«7> 


13/2,1/2)12,1) 


\n) 


|1,1>|0, 1/2) 




13/2,-1/2)12,2) 


\Te) 


|1, 1)11,-1/2) 



jiS*, Afs) |mij , m.i^) where S and Ms are the quantum 
numbers of the total electronic spin and its projection 
along the magnetic field. Furthermore, and nii^ are 
the nuclear spin projections of the individual atoms. In 
the MB the threshold energy E^{B) is equal to the Zee- 
man energy of the two atoms. Depending on the value 
of 5, the channel potentials correspond to the singlet 
(5 = 0) or triplet {S = 1) potential, i.e., V^{R) = Vs{R). 
While in the AB the coupling W^i^ is strong for small 
internuclear distances, in the MB the channels are only 
coupled by the relatively weak hyperfine interaction. The 
coupling W^,^ = (C'lV^f + V^^IC) is, on the other hand, 
present for all internuclear distances, which makes it im- 
possible to define open and closed channels in the MB. 

Depending on the distance between the two particles 
the set of interacting states is preferably considered in ci- 
ther of the two bases [3, . The AB of asymptotically 
uncoupled states is convenient for the description of the 
long-range part of the wave function. The MB is suitable 
for the short-range part where the exchange interaction 
leads to a strong coupling in the AB. While inappropriate 
for large distances, the MB is the natural choice to study 
molecular processes, such as the association of molecules, 
which take place when the atoms arc close to each other. 
Presently for ^Li-^^Rb the transition from the descrip- 
tion in the AB to the MB is appropriate at a distance 
i?sh ~ 20 ag (oq is the Bohr radius) where the exchange 
interaction is equal to the hyperfine interaction, i.e., 
where A£;M(^Li) + AShf(^'^Rb) = Joi?° -^e""-^, with 
AE^hfC^Li) = 228.2MHz and I^ii^'^Rh) = 6834.7MHz 
being the hyperfine splittings [39j . 



B. Computational details 



class of atomic and molecular collision problems. The 
essential idea is to divide the physical space into two 
or possibly more regions. In each region the stationary 
Schrodinger equation may be solved using techniques de- 
signed to be optimal to describe the important physical 
properties of that region. The solutions and their deriva- 
tives are then matched at the boundaries. The transition 
from AB to MB is carried out by a unitary basis trans- 
formation. 

The wave function \I' in Eq. ([3]) must obey appropri- 
ate boundary conditions in order to reduce the number of 
the independent solutions of the set of equations in ^ to 
one. The condition V'q(O) = ensures that the full wave 
function does not diverge at i? = 0. Another demand is 
that functions of the closed channels V'x(^) must vanish 
&t R oo. The implementation of these boundary con- 
ditions allows to solve Eqs. (|4]) leaving one free parameter 
in the solution, e. g., the normalization of the open chan- 
nel. We chose to scale the open channel function to the 
sin-normalized form 

V'ai(i?)|/?^oo =sin(fc-i? + (5), (9) 

with k = \J2[lE. The phase shift (5 is a result of the 
interaction and is connected via 

tan((5) = — fc • Csc , (10) 

to the s-wave scattering length Osc- In order to normal- 
ize the incoming channel function its asymptotic form is 
matched using Eq. ([9|). The value of Osc is automati- 
cally determined by the matching procedure. As will be- 
come evident in Sec lIIDI a variation of the magnetic field 
around a resonance leads to a transition of the phase 
through 7r/2 and thereby drastically changes the value 
of Osc- There arc different types of normalization, e.g., 
the energy or momentum ones. For calculating observ- 
ables like absolute transition rates the norm plays a role. 
However, general conclusions of the present work do not 
depend on the choice for the normalization. 

The kinetic energy E of two atoms when they are far 
apart is set to the arbitrarily chosen small value of 50 Hz. 
Since this energy is very small, the collisions are limited 
to the s-wave type only. The choice of a small but finite 
energy is justified because under ultracold conditions two 
particles collide with a low but non-zero energy. Further- 
more, the non-zero energy helps to avoid non-physical 
numerical artifacts in the definition of the phase i5. 



Since for the present case of ^Li-^^Rb the channel with 
the lowest threshold energy \a\) = |l/2, 1/2)|1, 1) is con- 
sidered as the open entrance channel, only channels with 
the total angular momentum Mp = 3/2 are coupled dur- 
ing the collision. All eight coupled atomic and molecular 
basis states are given in Tab. HI 

The system of eight coupled equations is numeri cally 
solved in the AB employing the i?- matrix method [34| . 
This method is a general ah initio approach to a wide 



C. Multi-channel results 

The system of ^Li-^^Rb features for a collision energy 
i? = 50Hz two s-wave resonances in the range of i? < 
1500 G, a broad one a,t B = 1066, 917 G and a narrow 
one at B = 1282.576 G (see Fig. [1]). While the narrow 
resonance is also examined, this paper focuses on MC 
solutions around the broad resonance. This resonance 
has been also observed experimentally [46| and is well 
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FIG. 1: (Color online) Scattering length function of the 

external magnetic field value B for ^Li-*^Rb scattering at i? = 
50 Hz. A broad and a narrow MFR are visible at Bq = 1066, 917 G 
and Bq = 1282.576 G. The horizontal line marks the background 
scattering length a^g = — 17.8ao of the left resonance. 
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FIG. 2: (Color online) The channel functions ip^{R) for the ^Li- 
^''Rb collision in an off-resonant field B = 1000 G (a) and a field 
B = 1066.9 G close to the resonance (b). The atomic labels (see 
Tab. |I]| are indicated in (a). The insets focus on a region of small 
internuclear distance. 
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FIG. 3: (Color online) The channel functions ?/)j(i?). The same as 
Fig. [2] but in MB. The molecular labels are indicated in (a). 



reproduced by the MC calculations. Moreover, processes 
like, e.g., PA are more efficient for a broad resonance 
because three-body losses can be minimized in this case. 

Figures [2] and [3] present the channel functions of the 
MC calculations in AB and MB, respectively, for a colli- 
sion of ^Li-^^Rb at two different magnetic field strengths 
B. Figures, ^b.) and [3l^a) show the case of a far off- 
resonant field of 5=1000 G which results in a small 
scattering length of only Usc = — 14.9ao. Figures [^l^b) 
and[3]^b) are taken close to the resonance at i?=1066.9G 
with a scattering length of age = —65 450 oq. This 
large value is arbitrarily chosen for the present study. 
It is already a good representation of the resonant case 
Use = -oo. 

The change of the long-range behavior between two 
scattering situations with small and large Osc can be more 
clearly analyzed in the AB where all but one channel 
are closed, i.e., decay for large internuclear separations 
(see Fig. [5]). As is evident from Figs. [D^a) and (b) the 
open-channel wave function ipai changes the slope result- 
ing in a plateau when changing Osc from small to large. 
Furthermore, the resonant open-channel function has a 
much larger amplitude within the considered range of in- 
teratomic distances than the off-rcsonant one. This large 
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difference in amplitudes (about four orders of magnitude) 
sustains in the region of small intcrnuclear distances (see 
insets of Figs. Elja) andEIb)). 

At small intcrnuclear distances above i? w 7 ag the 
channel functions in the AB show quite irregular behav- 
iors (see insets of Fig. , which is a result of the large 
coupling proportional to the exchange energy Vc:k{R). In 
the MB the coupling between the channels is induced by 
the hyperfine interaction that is much smaller. Hence, 
the channel functions show a clear behavior of pure sin- 
glet and triplet wave functions for small internuclear dis- 
tances (see insets of Fig. For distances R < 7 qq the 
triplet components vanish due to their higher exchange 
energy. Accordingly, also in the AB the channel functions 
are similar to pure singlet wave functions ai R < 7 uq (see 
insets of Fig. [2]). All channel functions in the MB con- 
tribute correspondingly to the decomposition of the open 
channel ipai into states of the MB. Therefore, at large 
internuclear distances they look similar to ipai ■ It is im- 
portant to note that at small distances the closed channel 
functions have non-zero amplitudes even in the B-field- 
free case; they are slightly excited during the collision and 
possess a background contribution to the scattering pro- 
cess. Therefore, the two-body collision is a multi-channel 
process even in field-free space. 

Due to the resonant coupling at i?=1066.9G, the ad- 
mixture of the closed channels increases about four or- 
ders of magnitude. This is well described by the TC ap- 
proximation [2^, [2^ where the admixture of the closed 
channel and the long-range behavior of the open chan- 
nel show a similar dependence on the scattering length 
(see Sec. Ill Pp . In contrast to the TC approximation 
where one assumes that a bound state composed of a su- 
perposition of all closed channels is simply scaled at the 
resonance, the relative amplitudes change in reality be- 
tween the resonant and off-resonant cases. On the other 
hand, the functional form of all closed channels indeed 
stays constant (compare, e.g., channel |a4) in Figs.[2I^a) 
andlUb)). Altogether, this gives hope to be able to re- 
produce the change of the amplitude of both the open 
channel and the closed channels at small internuclear dis- 
tances around an MFR with just one SC wave function. 



D. Two-channel approximation 

The TC approximation is very successfully used to de- 
scribe resonance phenomena in MC problems [20l. [25l. [2l| . 
It is briefly introduced in order to understand to what 
extent SC approaches can mimic MC systems. A more 
rigorous introduction may be found in (33l . IZtI . |48| . 

Within the TC approximation one projects the MC 
Hilbert space onto two subspaces, the one of the closed 
channels (with projection operator Q) and the one of 
the open channel (with projection operator P). The full 
wave function is thus written as j^*) = (P + Q)|^') = 
l*p) + I^q)- An MFR occurs, if the energy E of the 
system is close to the eigenenergy Eq{B) of a bound 



state of the closed-channel subspace. In the one-pole 
approximation one effectively assumes that the closed- 
channel wave function is simply a multiple A of the bound 
state l^b), i.e., |^'q) = This approximation 

yields the closed- channel admixture [sl] 



(11) 



A = -C\l — sinJrc 



where C is a normalization constant. The long-range 
behavior of the open channel is given as 



*p(^)lfl^oo = ^\/ ^ sin(fci? + <5bg + Jrcs) . (12) 



If the wave function is sin normalized, then C = 
Another popular choice is the energy normalization with 
C = 1. However, the presence of an external trap can 
also induce a dependence of the normalization on the 
long-range behavior of the open channel parameterized 
by flsc, such that in general C ~ C(asc). 

The total phase shift S — (5bg + ^ros results from the 
background phase shift 6hg of the open channel without 
coupling to the closed channels and from a contribution 
(5ros due to the resonant coupling to the bound state. 
Via Eq. pO)) the total phase shift is connected to the 
scattering length Qsc- The TC approximation yields for 
k the well known relation 14911 



abg 1 



AB 

B- Br 



(13) 



between scattering length and magnetic field strength, 
where Obg = — tan 6hg/k is the background scattering 
length, AB is the width of the resonance, and Bq its 
position. 

We note that independently of the normalization func- 
tion C(asc) both the admixture of the closed channel 
A and the long range open-channel solution (fT2|) show 
for small energy, not too large internuclear distances 
(i.e., kR ^ S), and small backgroimd phase shifts (i.e., 
S w iJrcs) a similar dependence on the scattering length 

^sc ■ 

Usually, for small energy E the background phase shift 
<^bg = — arctan(A:abg) is necessarily also small. Since a 
scaling of the open-channel wave function in the long 
range is more or less directly continued to shorter dis- 
tances, the proportionality between A and ^p{R) holds 
approximately also for smaller R. Therefore, looking at 
molecular processes taking place at small internuclear 
distances, the enhancement of the closed- channel contri- 
bution is already reproduced by the open channel. This 
paves the way to an SC description which will now be 
discussed. 
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III. SINGLE-CHANNEL APPROACHES 
A. Variations of the single-channel Hamiltonian 

In order to reflect the molecular behavior at small dis- 
tances, we will seek to base the SC approximations on 
pure singlet or triplet interaction potentials. This en- 
sures that the nodal structure of the resulting SC wave 
function is similar to the relevant singlet or triplet com- 
ponents of the MC system. The final aim is to mimic 
in parallel the long-range behavior of the open channel 
and the variation of the amplitude of singlet or triplet 
components in the vicinity of an MFR. 

In an SC approach the interaction strength can be arti- 
ficially varied by a controlled manipulation of the Hamil- 
tonian 



1 



V{R). 



(14) 



Subject to modification arc the inter-atomic potential 
V{R) and the reduced mass /i of the system. The mod- 
ifications can lead to a shift of the energy of the least 
bound state relative to the potential threshold. When 
lifted above the threshold, the bound state turns into a 
virtual state [i^, \E^. A large scattering length of the 
solution of the SC Schrodinger equation with Hamilto- 
nian can be elegantly explained by a resonance of 
the scattering state with either a real bound state or a 
virtual state close to the threshold [H, HO] ■ Within an 
SC approach the energy of a bound or virtual state is 
changed in order to induce a variation of the scatter- 
ing length. In this respect SC approaches show striking 
similarities to MFRs where the energy of a bound state 
in the closed-channel subspace is moved by changing its 
Zeeman energy by an external magnetic field. 

As argued before, the SC wave functions should be ei- 
ther of singlet or triplet character for small internuclear 
distances. Wc reduce our considerations for Li-Rb to the 
singlet case and chose as initial potential the one for the 
electronic ground state, i.e., V{R) = Vxi^+{R)- 
This potential is varied by a controlled manipulation of 
the strong-repulsive inner wall |lll|. the long-range van 
der Waals attraction Vvdw {R) [23] , and a novel Gaussian 
perturbation around the transition point Rgh between the 
molecular and the atomic description of the system (in- 
troduced in Sec. Ill A"2)) . These procedures will be called 
s variation, Cg variation, and G variation, respectively. 

The potential variations are induced by replacing V{R) 

by 



\viR) 



-Re --Re 



V"^{R) = F(i?) +(5GExp 



) R<Re 
R > i?e 



R - Rg 



or 



V'^^{R) = V{R) + ^-f{R) 



(15) 



(16) 



(17) 



where Re = 6.5 ao is the equilibrium distance and Rc = 
4.6 oq is the crossing point of the Vxi-s+{R) with the 
threshold. The width in the G variation is chosen as 
CT = 2 ao and its position as Vq = Rsh + f. The smooth 
variation of the long-range region of the potential in the 
Cg variation is achieved by the gradual stepping function 



f{R) 



1 



(18) 



where 7 — ln(999) w 6.9 ensures that f{R) rises from 
0.001 to 0.999 in the region Rq - A < R < Ro + A. For 
the present study the parameters of the tuning function 
are chosen as A = 6 and Rq = IGuq. The three potential 
variations are depicted in Fig. 51 
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FIG. 4: (Color online) Original XiS+ potential V{R) (blue solid) 
with applied s variation (black solid), Cq variation (dashes) and G 
variation (dots). The variation parameters are s = 0.03 ao, SCa = 
Cfj/2 and 5G = ^(-Rijij) (see Tab. HTlfor the adopted parameters). 
The insets show some relevant ranges of R on an enlarged scale. 

An alternative way to tune Qsc is offered by the /i vari- 
ation within which one changes the reduced mass of the 
system by /z ^ fi — [13]. This alters the kinetic en- 
ergy operator and can modify the energy of the least 
bound state like potential variations. In contrast to the 
presented potential variations which act on cither the 
short-range, mid-range or long-range part of the poten- 
tial, the mass variation influences the Schrodinger equa- 
tion at any distance. It is very similar to a scaling of the 
potential by V{R) -fV{R) [H]. The only difference is 
an additional change of the energy-momentum relation 
E{k) = fc^/(2/x) which can, e.g., slightly influence the 
normalization of the wave function. 

One can think of several other approaches to vary the 
SC Hamiltonian. For example, one can vary the strength 
of the exchange energy Jq, its decay parameter a or the 
van der Waals parameters Cs, Cio |27l . [5l| . The current 
approaches are chosen to comprise variations that act on 
the short rage of interatomic distance (s variation), on 
an intermediate range (G variation), on the long range 
(Cg variation), and on the full range {fj, variation). 

No matter which SC approach is flnally chosen, a map- 
ping between the MC system and an appropriate SC 



TABLE II: Values of the parameters for s,Ce,G,iJ, varia- 
tions at small and large scattering lengths resulting from 
a magnetic field far away, close and right at an MFR. An 
infinitesimally small change of the field right at the reso- 
nance {B = 1066.92 ao) switches the interaction regime from 
the infinitely attractive (ttsc = —oo) to infinitely repulsive 

(flsc = +00). 



B /Gauss 








<5G/|V(i?,h)| 


5fi/fi 


1000.00 


-14.93 


-0.00947 


-0.0281 


-0.1310 


-0.00235 


1066.90 


-65450 


-0.04142 


-0.1282 


-0.6296 


-0.01046 




—00 


-0.04145 


-0.1283 


-0.6302 


-0.01047 


1066.92 


+ OQ 


0.13065 


0.3357 


-4.3856 


0.02984 



Hamiltonian is straightforward. Knowing the parame- 
ters AB and Bq in Eq. (fT^ for an MFR cither from 
experimental data or a coupled MC calculation one can 
connect each value of the magnetic field B to a scattering 
length Gsc and a corresponding value of the SC variation 
parameter that induces the same value of the scattering 
length. Clearly, this additional information is required, 
i. e., the SC model has no predictive power by itself. 

Typical values of the four variation parameters as 
they will be used in the following are given in Tab. [Tll 
The wave functions resulting from the different varia- 
tion methods are denoted (j)^'{R) where v G {s, G, Cg, ^} 
stands for the applied v variation. 



B. Multi-channel vs single-channel 

In the following, the wave functions of the SC and MC 
approaches are compared. As discussed in Sec. Ill A 2l the 
appropriate choice of the MC basis depends on the inter- 
atomic distance. While for large interatomic distances 
{R > Rsh) the description in the AB is adequate, their 
basis states arc strongly coupled for shorter distances. 
Here, the MB describes the physical properties far bet- 
ter. Due to the weak hyperfine coupling the states in the 
MB keep to a good degree of accuracy the structure of 
an uncoupled singlet or triplet state, respectively. Close 
to an MFR solely amplitudes for some of the states are 
heavily increased. Figure [5] shows for example a com- 
parison of the singlet state \Si) close to an MFR with 
the same state far away from the resonance and with 
the other singlet state |52) again close to the resonance. 
Clearly, for R < Rgh they differ only by a constant pref- 
actor. This is important, as it allows to describe the 
short-range behavior of the MC wave function by either 
a pure singlet or triplet channel function depending on 
the physical process which is to be described. For ex- 
ample, only singlet components contribute to the DPA 
process for the transition into the absolute vibrational 
ground state (as it will be discussed in Sec. lIVp . hence, 
the triplet components may be omitted. 

In the following, the aim of the SC approach is to 
mimic the behavior of the MC singlet components for 
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FIG. 5: (Color online) The channel functions of the singlet state 
\Si) close to the resonance (age = —65 450 ao) and away from the 
resonance (osc = — 14.9ao) are depicted together with the channel 
function of the close to resonant IS2) state. All three functions 
differ for R < 30ao only by a constant prefactor. 



R < Rsh by a controlled variation of the SC Hamilton 
operator (fH|) with singlet potential Vxi^+{R)- With 
the help of the s,Ce,G and fj, variations presented in 
Sec. nil Al the SC wave function is adjusted to match 
the asymptotic behavior (i.e., the scattering length Csc) 
of the open channel for a given external magnetic field 
B. The cases of an off-resonant magnetic field {B = 
1000.0 G) and one close to a resonance {B = 1066.9 G) 
are considered. The corresponding scattering lengths are 
ttsc = — 14.9ao and age — —65 450ao (see Sec. Ill C] and 

Figs, o ED. 

Figure [5] shows a comparison of the SC wave functions 
(j)'" (R) with the channel functions of the dominant singlet 
channel \Si) and the open channel \ai) for the full range 
of short and long interatomic distances. Figure [B] allows 
to examine how the different variational methods are able 
to reflect both the behavior of the singlet components for 
distances R < Rsh and the one of the open channel for 

R > Rsh- 

Generally, any SC approach has to induce a shift 
of the phase S in order to tune the scattering length 
asc = — tan((5)/fc. The difference 6 — (5ini from the 
phase of the unperturbed system i5ini is accumulated 
where the variation of the SC Hamiltonian takes place. 
Since the scattering length of the original singlet poten- 
tial VxiE+(-R) is with Og"' = 2.3 ao relatively close to 
Qsc = — 14.9ao7 hardly any phase shift has to be ac- 
quired {S — (5ini = 8.6 • 10~^7r) to match the open channel 
for the off-resonant magnetic field B = 1000 G. Accord- 
ingly, the nodal structure of the MC singlet component 
is very well matched in Fig.[6Ua). The situation changes 
close to the resonance where the large scattering length 
Use = —65 450 ao requires a phase shift of S — Sini ~ 0.227r 
(Fig. El^b)). This is about half way to the resonant phase 
shift 7r/2. 

For the s variation the total phase shift to match the 
open channel is acquired for distances R < 6.5 ao. Ac- 
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FIG. 6: (Color online) Comparison of the SC wave functions <p^'{R) 
with the MC functions of the singlet state (scaled) and the 

open channel |ai). The SC potentials are varied to match the 
asymptotic behavior of the MC channel functions of the open chan- 
nel, (a) Off-resonant case with asc = — 14.9ao {B = 1000 G). (b) 
Resonant case with Qbc = —65 450 ao {B = 1066.9 G). The accord- 
ing values of the s, SCe, SG, and Sfi parameter are given in Tab. |lT] 
The smaller plots focus respectively on a region of small internu- 
clear distance (left) and a region R ij^ij (right). 



cordingly, the nodal structure between the channel 
function and the SC (p^iR) wave function is shifted for 
R > 6.5 oo (see upper left plot in Fig. [B]Jb)). Contrar- 
ily, both the Cg variation and the G variation induce a 
phase shift for distances R larger than Rq — a = 20 cq 
and i?o ~ IGflQ. Thus, for smaller internuclear distances, 
the SC wave functions 4)'~^{R) and <t)'~^^{R) coincide with 
the l^i) channel function. Finally, since the ^ variation 



acts on any internuclear distance, the phase difference is 
gradually accumulated for (j)'^{R). 



Depending on the range of variation for the SC ap- 
proaches, also the matching to the open channel of AB 
differs. The s variation matches the open channel al- 
ready closely before Rsh- Surprisingly, also the fi vari- 
ation shows a reasonable match already before Rsh, al- 
though it acts also for larger distances by changing at 
least the dispersion relation E{k). This effect may, how- 
ever, not be visible, since kR <C 1 in the plotted region. 
The Cg variation changes the long-range behavior of the 
interaction potential. Correspondingly, the wave func- 
tion shows a clear difference to the open channel even up 
to i? = 100 ao. The Gaussian perturbation of the G vari- 
ation acts only around Rsh_. This results in the favorable 
situation that both the |S'i) channel for R < 20 and 
the open channel for i? > 24ao are matched by the SC 
wave function. 



Since the nodal structure among different singlet and 
different triplet channels coincides for R < Rsh the pre- 
sented results are generalizable to any singlet or triplet 
state. Thus, SC approaches are generally able to repro- 
duce the asymptotic behavior of the open channel of the 
MC wave function in the presence of an MFR while also 
reflecting certain aspects of singlet or triplet components 
for small internuclear distances. Depending on the re- 
gion of the variation of the SC Hamiltonian, the nodal 
structure of any channel function in the MB can be re- 
produced for R < Rsh- The most flexible SC approach is 
the G variation which is able to smoothly switch between 
the accurate description of a MB channel and the open 
channel. Furthermore, it offers the advantage, that one 
can define the transition point (here R = Rsh) at will, 
such that also for slightly larger distances MB channel 
functions can be emulated. 



An aspect of the MB channels which cannot be re- 
flected by the present approaches is their absolute am- 
plitude. Since the amplitudes at small internuclear dis- 
tances of the different channels change drastically in the 
presence of an MFR, they have a large impact on molecu- 
lar processes such as the association of molecules utilizing 
MFRs. In the next section the exemplary case of a di- 
rect dumping of the scattering state to the vibrational 
ground state of the AT^E^ is considered. The transition 
rate depends strongly on the behavior of the amplitude 
of the dominant singlet state \Si) which was considered 
in this section. It will be shown that although the ab- 
solute amplitude of this state is not reproduced by any 
SC approach, the relative enhancement of the transition 
rate at magnetic fields close to a resonance can be well 
refiected. 
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IV. PHOTOASSOCIATION OF "LI-^'^RB TO 
THE ABSOLUTE VIBRATIONAL GROUND 
STATE 



Ultracold polar molecules are of great interest for many 
applications in quantum information processin g [ssl . |36| . 
the exploration of lattices of dipolar molecules 37 1 , pre- 
cision measurement of fundamental constants 5^, and 
ultracold chemical reactions [sl, [U- Since standard 
cooling technics developed for atoms are not suitable 
for molecules due to their complex level structure, ultra- 
cold molecules may alternatively be associated directly 
from ultracold atoms. As was already mentioned in 
the introduction the starting point to create ultracold 
molecules in their vibrational ground state are often Fesh- 
bach molecules formed by a sweep of the magnetic field 
around an MFR in a high- lying vibrational level [20| . 
These loosely bound molecules are usually transferred 
by complex PA schemes via intermediate excited states 
to the desired vibrational grou nd state [sl, [s^. Espe- 
cially STIRAP [13, S showed to be successful 
in efficiently creating ultracold ground state molecules. 
However, Feshbach molecules possess a relatively short 
life time such that a Feshbach optimized transition di- 
rectly at the resonance can be favorable ^ . 
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FIG. 7: (Color online) Sketch of the resonant SC wave functions 
obtained via G variation and respective BO potentials relevant for 
the DPA transition to either the singlet ground state (thin, solid) 
or the triplet ground state (thin, dashed) of the respective X^S+ 
(thick, blue) and a^S+ (thick, red) potentials. For better visibil- 
ity, the potentials and wave functions arc shifted along the y-axis. 
In reality, singlet and triplet potentials have the same threshold 
energy. 



Calculation of transition rates 



For all schemes that make advantage of the resonant 
coupling to a molecular bound state at an MFR [sl. [lOl. [2^. 
IHj El, E^i HEl , the increase of the amplitude for the rel- 
evant channels as the scattering length grows is of great 
importance to enhance the molecule creation. Although 
in the last section it was shown that the absolute am- 
plitude of the MB channels is not reproduced by the SC 
approaches, the TC approximation gives hope that the 
relative enhancement can still be recovered. In Sec. Ill Dl 
it was discussed that both the admixture of the closed- 
channel bound state and the open-channel function scale 
similarly with the scattering length. One can therefore 
expect to be able to combine this collective relative en- 
hancement into one channel. 



In the following, the Feshbach optimized DPA 
(FOPA) [6l| to the absolute vibrational ground state of 
^Li-^^Rb in the electronic X^E"*" state is considered to 
examine the applicability of SC approaches to study pro- 
cesses of molecule creation. We consider this case since it 
has an interest on its own for the creation of bound ultra- 
cold molecules. Furthermore, the transition rate to the 
absolute vibrational ground state depends on the scat- 
tering wave function at very small internuclear distances 
(see Fig. [7|) . We also examined the transition to the vi- 
brational ground state of the electronic triplet state a^'S'^ 
which is situated at slightly larger interatomic distances. 
Since we found no essential differences to the singlet case, 
we focus on presenting only its results in this work. 



Given the solution of the MC problem ^{R) = 
^''r^'^ 10 the MB for a given magnetic field B 
the free-bound FOPA transition rate ^i{B) to the final 

molecular state *f(i?) = ^^^y/(e, $)|Cf) with vibra- 
tional quantum number v and rotational quantum num- 
ber J within the dipole approximation is proportional to 
the squared dipole transition moment [62l | 



'MC 



■^^{R)D{R)'ilj^,{R)dR 



(19) 



Here, D{R) is the electronic dipole moment. Within the 
dipole approximation only transitions from the s-wave 
scattering function to a final state with J = 1 arc allowed. 
Due to the orthogonality of the MB, only one molecular 
channel has to be taken into account in Eq^ (fTT ^ 
The TC approximation predicts a rate [32 



Itc{B) = \C-Cf-\sm{S,,,{B)-6o)f 



(20) 



where the constants C and Sq, explicitely given in [33[, do 
not depend on the magnetic field within the TC approxi- 
mation. The phase shift 5o is usually small [33| and thus 
the minimum lies close to a vanishing resonant phase shift 
(5rcs = 0, i.e., close to the background scattering length 
Obg- We determine Obg by a fit of asc{B) to Eq. ([13]) 
which yields Obg = —17. Sag. From S — (5res + ^bg and 
Eq. ([Tn|) one can then directly determine 5yos{B). The 
behavior of Eq. ^U\\ accurately reflects the one of a MC 
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system for well separated resonances [SSj] . We use it here 
to determine the maximal MC transition rate. 

The transition rate to the final state within an SC 
approach is simply proportional to 



^,{R)D{R)cf>''{R) 



where v, as before, denotes the variational method which 
for the present analysis induces the scattering length age 
equal to the one of the MC system for a given i?-ficld 
value. 

D(R) is again the electronic dipole transition mo- 
ment. For the purpose of the present study we reduce 
our considerations to the linear approximation D(R) = 
Dq + Di- R. The SC scattering wave function is orthogo- 
nal to the different vibrational bound states. In the MC 
case only the weak hyperfinc coupling in the MB causes a 
very slight non-orthogonality. The influence of Dq can be 
therefore safely ignored. Calculations with higher-order 
expansions showed that the exact functional behavior of 
D{R) (obtainable from [6^) does hardly influence the 
relative enhancement of the transition rate. Thus, the 
use of D{R) = Di ■ R does not restrict generality. It 
is important to note that Eqs. P^ - (PT|) are only valid 
within the dipole approximation. It is supposed to be 
applicable, if the wavelength of the associating photon 
is much larger than the spatial extension of the atomic 
or molecular system. The shortest PA laser wavelength 
corresponds to the transition to the lowest vibrational 
state. Although the spatial extention of the initial state 
is infinite, the integrals for dipole transition moments is 
finite, as it contains a finite wave function of the bound 
vibrational state as a factor. Therefore the dipole ap- 
proximation is valid. 



B. Comparison of transition rates 

A change of Ogc leads to an increase or decrease of F". 
In order to quantify the magnitude of this change, an en- 
hancement or suppression factor may be introduced 



(22) 



It describes the relative enhancement [g"" > 1] or sup- 
pression [g" < 1] of the DPA rate at a given Osc vs. a ref- 
erence scattering length alf , for a specific final state i;. 
Although it may appear to be most natural to choose 
a™^ = 0, a large non-zero value offers some advantages. 
In this case, /"(a^^^) is not too small and large numerical 
errors are avoided. 

Figure [5] shows a comparison of the SC transition rate 
for the different variational approaches with the correct 
MC result. In the calculation of the MC transition rates, 
we assume a measurement in which the nuclear spins are 
not resolved. This corresponds in practice to the case in 



which the transition rates from the \Si) and IS'2) channel 
are summed. All rates are normalized to their respective 
maximum value {al^^ — 00). Note, however, that the 
different absolute dipole transition moments disagree by 
some orders of magnitude. 
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FIG. 8: (Color online) Comparison of MC and SC results for the 
transition rate to the absolute vibrational ground state relative to 
the respective maximal transition rate as a function of the scat- 
tering length (a) and the inverse scattering length (b). The MC 
results arc fitted according to the TC approximation (Eq. (I20I I). 

A fit of the MC result by the TC estimate with only two 
free parameters C and 60 reveals that the simple depen- 
dence of the transition rate given by Eq. (pO)) describes 
the transition process of the MC system correctly. 

All SC approaches agree with the MC result for large 
scattering lengths in the proximity of the resonance (see 
Fig. EUb)). For small scattering lengths where the tran- 
sition rate is already suppressed by more than four or- 
ders of magnitude deviations from the MC result appear. 
The differences mainly originate from a shift of the min- 
imal transition rates of the SC approaches compared to 
the MC result. In the MC case the minimum lies at 
Use ~ — 21.1ao close to the background scattering length 
Cbg = — 17.8ao in accordance with the TC approxima- 
tion. The minima of the SC approaches tend to be sit- 
uated on the positive side around Qsc ~ 50aQ. This is, 
however, not a general trend, since we observed for other 
transitions also minimal SC transition rates at negative 
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scattering lengths. The location of the minimum depends 
on a system under investigation and on the applied SC 
variation. 

Figure [5]^a) features two kinks of the transition rate 
at Osc = 40 ao for the s and fi variations. This can 
be explained by the shift of the nodes of the SC wave 
functions which takes place at the equilibrium distance 
of the bound molecule and therefore influences the PA 
rate. Since the variation parameters are tuned around 
their resonance value, with increasing distance from the 
resonance both left and right of it, eventually the same 
scattering length is induced (see Fig. [S]). However, the 
nodal structure of (f>''{R) and 4'^{R) for short ranges can 
differ, leading to different transition rates. This does not 
occur for the Cq and G variations that act far beyond the 
equilibrium distance. Note, however, the scale at which 
the kink is visible. Its effect on the rate is minute. 

Analogous examinations were also done for the other 
MFR of ^Li-**^Rb at B = 1282.58 G. Although this reso- 
nance is two orders of magnitude narrower than the one 
considered before and the amplitudes of the channels are 
different, no significant differences for the relative rates 
were observed. The generality of our considerations is 
also supported by calculations of the dumping rate to 
the vibrational ground state of the triplet configuration 
o?Yj'^ . In all cases the SC approaches showed a compa- 
rable ability to reflect results of the MC system. 

It is also interesting to note that results of the anal- 
ysis show that neither the details of the interatomic nor 
magnetic-field interactions are relevant for the calcula- 
tion of the relative rate. A simple SC model turns out 
to be adequate to calculate the relative enhancement of 
the PA process. Furthermore, in view of the important 
question of how to optimize the efficiency of DPA, Fig. [H] 
reveals once more that the use of a large absolute value 
of the scattering length is favorable. 



C. Number of bound states 

As already mentioned, SC resonances are evoked by ar- 
tificially shifting the least bound state or a virtual state 
across the threshold. By turning a bound state into a vir- 
tual state or vice versa, the total number of bound states 
TVb changes necessarily. This can be avoided by stopping 
the variation just before the bound or virtual states reach 
the threshold. Nevertheless it is possible to achieve any 
scattering length by moving between two different SC 
resonances. This is illustrated by the example of the /i 
variation in Fig. [D^a) where three resonant branches of 
the asc{Sfi) curve are depicted. As discussed in [l^ the 
question arises, whether it is preferable to keep iVb con- 
stant or to change the variation parameter across a SC 
resonance as was done so far in this work. 

In Figs. [TOf a) and (b) the relative transition rate is 
depicted as a function of the phase shift 6. In compar- 
ison to a l/osc-plot (Fig. [5Jb)), this allows an enlarged 
view on the region of resonance where the phase S sud- 
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FIG. 9: (Color online) (a) Scattering length Qgc clS SL function of the 
(5/i-parameter of the mass variation. By constraining the variation 
to the thick (red) branch, any scattering length is reached while 
keeping the number of bound states N^, constant. By constraining 
it to the dashed (black/red) branch any scattering length is reached 
while A''i3 changes, (b) Resonant SC functions [asc = oo) at two 
different (5/i-parameters 5fi fa — 0.01/j (dashed), 5fi O.OSfi (thick, 
red). In order to make the relevant phase and amplitude difference 
at small internuclear distances visible one of the wave function is 
multiplied by —1 in the inset. 



denly crosses tt/2. In Fig. fTUTal the SC variations are 
performed in the same way as in Fig. [H] around one SC 
resonance while changing A^b- This results in a perfect 
agreement with the MC result (the deviations at small 
relative rates as shown in Fig. [H{a) are not visible on 
a linear scale of the relative rate). Furthermore, large 
values of the scattering length can be obtained by slight 
modifications of the SC Hamiltonian. By fixing Ab, one 
has to stay on the same branch of the resonant curve 
Usciv). This modifies the SC Hamiltonian strongly and 
can lead to a sudden change of the relative rate by some 
30% as is shown in Fig.ITOJb). 

The reason for the sudden change of the wave function 
is twofold. In all cases the asymptotic behavior of the 
wave functions are the same at different resonant points 
corresponding the same age, but different Hamiltonians 
lead to a slightly different continuation of the wave func- 
tion towards smaller distances. If the variation takes 
place at ranges larger than the equilibrium distance Re 



13 



(a) 



1.0 








0.8 






■ 


0.6 






■ 


So 




/ • • • MC result \ 




0.4 




/ ^— Fit of MC result \ 
/ s vsri^tion \ 


■ 


0.2 




/ Q variation \ 

/ .... vcii'icil'irtn \ 

// variation 


. 


0.0 













4 2 4 

Phase shift 6 


n 


(b) 









1 1 1 1 1 1 1 1 r 




—1 1 1 1 1 1 1 1 1— 




MC result \ 
Fit of MC result\\ 
s variation 
Cf, variation 

G variation V- 
fi variation 
J I 1 1 1 1 l_ 



Phase shift 6 



FIG. 10: (Color online) Comparison of MC and SC results for the 
transition rate to the absolute vibrational ground state relative 
to the respective maximal transition rate as a function of the phase 
shift 5. The different SC variation parameters are cither varied 
around one SC resonance (a) or between twro resonances staying on 
the same aaciSfj.) branch (b). The MC results arc fitted, according 
to the TC approximation (Eq. I I20I I). Note the sin^-like form of the 
functions. 



(Cg , G and fi variation) , the wave function around Re 
where it influences directly the transition rate can dif- 
fer shghtly in ampUtude. Secondly, if the variation takes 
place around Re (s and /i variation) the nodal structure 
of the SC wave functions differs for both resonant SC 
parameters, since the necessary phase shift is acquired 
in different ways. Both effects induce a "step" in the 
transition rate at S = 7r/2 as is visible in Fig. [TUT b). 
The nodal shift can in principle change the transition 
rate more strongly than in the present case. Figure [5I^b) 
compares the two wave functions of the fi variation at dif- 
ferent resonant S/j, parameters. One can observe around 
R^ Re both effects just described: the change of ampli- 
tude and the change of the nodal structure for different 
resonant variation parameters. 

To conclude, in order to calculate relative PA rates it 
should be in most cases preferable not to keep the num- 
ber of bound states fixed and to avoid a sudden change 
of the SC scattering wave function while going over the 



resonance. The drawback is of course a sudden change 
of the wave function for small scattering lengths. But 
here the SC approaches show in any way differences to 
the MC result, such as a shift of the minimal transition 
rate along the (Fig. IDJa)). Noteworthy, for the 

energy spectrum analysis, as it was done, e. g., in [ll| for 
two atoms in an OL, it is more convenient to stay on the 
same SC resonant branch. The alternative variation with 
non-constant A^b does not influence the resulting energy 
spectrum. However, the disadvantage is that the num- 
bering of the discrete levels should be changed across a 
SC resonance. 



V. CONCLUSION 

We presented single- channel approaches that were able 
to reproduce both the long-range behavior of the open 
channel as well as the nodal structure and relative en- 
hancement of any singlet or triplet state of a multi- 
channel system in the presence of a magnetic Fesh- 
bach resonance. However, single-channel variations in- 
duce a shift of the nodal structure not present in the 
multi-channel solution. Furthermore, the overall ampli- 
tude of the wave function stemming from the asymp- 
totical behavior can be slightly modulated by long and 
intermediate-range variations. The G variation, intro- 
duced in this work, showed to reproduce the correspond- 
ing multi-channel components at short and long inter- 
atomic distances most accurately. 

As was demonstrated for the exemplary case of ^Li- 
^^Rb scattering single-channel wave functions can be 
used to study processes of molecule formation. We exam- 
ined the specific process of a direct one-photon photoas- 
sociation to the absolute vibrational ground state of ^Li- 
^^Rb and proved the applicability of the single-channel 
approaches to model this process. The effects of the 
nodal shift and the modulation of the amplitude lead to a 
discontinuity in the transition rate for either small scat- 
tering lengths, if varying the single-channel Hamiltonian 
over a resonance, or at large scattering lengths, if keeping 
the number of bound states constant. As was discussed, a 
variation around a resonance of the single-channel Hamil- 
tonian is preferable, since the point of minimal transition 
at small scattering lengths deviates in any way between 
multi-channel and single-channel results. These devia- 
tions appear, however, at scattering lengths where the 
transition rate is negligible compared to the one at reso- 
nance. 

The general applicability of single-channel approaches 
was based on the two-channel approximation which re- 
veals that the scaling of the open-channel wave function 
and the admixture of closed channels depends on the 
scattering length in a similar way. Additionally, by the 
help of this approximation one is able to reproduce ex- 
actly the multi-channel transition rate by adjusting two 
free parameters, that combine all details of the transition 
process. 
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We can conclude that single-channel approaches are a 
suitable starting ground to study molecular processes in 
regimes were full multi-channel calculations are too labo- 
rious. This is, e. g., the case, if the scattering takes place 
in an external trapping potential like an optical lattice, 
that in general couples relative and center-of-mass mo- 
tions and spoils the spherical symmetry. In most cases 
the trapping potential does not directly influence the 
scattering wave function at short interatomic distances, 
but it induces an additional modulation of the amplitude 
as a function of the scattering length. The examination 
of effects due to these modulations are perfect candidates 
for the use of single-channel approaches. 

Since the nodal structure of either the singlet or the 
triplet components of the multi-channel wave function is 
reproduced by single-channel approximations, also more 
complicated photoassociation schemes, exciting a range 
of higher vibrational states, can be examined in the pres- 
ence of a trapping potential. Furthermore, single-channel 
approaches allow to treat three and many-body collisions 
with reasonable numerical efforts in the presence of a 
magnetic Feshbach resonance. 

Of course, the presented single-channel approaches 
have also clear restrictions. For example, one has to 
assume that the scattering energy and the background 
scattering length are sufficiently small. This condition 
can be spoiled for certain atomic systems and in deep ex- 
ternal trapping potentials with significantly large ground 
state energy. Another problem can be caused by the en- 
ergy dependence of the scattering length especially for 
narrow Feshbach resonances. This energy dependence 
is not reflected by the current approaches. Furthermore, 



the multi-channel wave function might behave differently 
compared to the single-channel one, if an energy varia- 
tion is induced by, e.g., ramping up an external trap. 
There exist single-channel approaches, which account for 
the energy dependence of the scattering length by a well- 
barrier pseudo-potential (63 |. However, like any pseudo 
potential it is unable to reflect the nodal structure of the 
scattering wave function at small intcrnuclear distances. 

Recently Deiglmayr et al. observed for ^Li-^^'^Cs (at 
5 = 0) the exceptional case of a strong deviation of 
molecular channel functions from pure singlet or triplet 
behavior at small intcrnuclear distances (26[ . Since spin- 
orbit coupling was neglected, they attributed this un- 
usual effect to strong hyperfine coupling but gave no 
reason, why this happens specifically for the considered 
system. It is certainly interesting to further investigate 
this effect which would limit the applicability of single- 
channel wave functions to predict, e. g., the relative tran- 
sition rates to different vibrational levels. 

Apart from this unusual behavior it should be possible 
from the theoretical considerations presented in this work 
to determine, whether and which single-channel approach 
is applicable for a specific system and molecular process. 
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